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ABSTEIACT 


Recently  an  efficient  technique  was  proposed  for 
performing  a  Fourier  transformation  on  integrands  that  are 
composed  of  aliased  factors.  In  this  report  the  proposed 
technique  is  utilized  to  perform  Fourier  transformations  on 
quantities  that  pertain  to  the  response  of  a  regularly  ribbed 
cylinder.  The  computations  relate  to  Fourier  transformations 
from  the  axial  wavenumber  domain  to  the  axial  spatial 
domain.  The  phenomena  of  aliasing  and  pass  and  stop  bands 
are  of  particular  interest.  It  is  argued  that  the 
complementarity  of  the  data,  when  presented  in  the  two 
Fourier  conjugate  domains,  may  be  a  useful  analytical  (and 
experimental)  tool. 
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1.  INTRODUCTION 


In  a  recent  report  the  authors  proposed  a  novel  computational  technique  to 
perform  a  certain  class  of  Fourier  transformations  [1].  The  class  pertains  to  Fourier 
transformations  in  which  the  integrands  are  completely  or  partially  composed  of  terms 
and  factors  that  are  aliased  with  respect  to  a  specific  scale  factor  in  the  domain  in  which 
the  integrand  is  expressed.  This  class  harbors  a  number  of  phenomena  in  structural 
acoustics.  A  few  of  these  phenomena  are  employed  in  this  report  to  exemplify  the  use  of 
the  double-sum  technique  for  computing  Fourier  transforms  [1].  In  particular,  the 
circumferential  modal  response  of  a  regularly  ribbed  cylinder  is  transformed  from  the 
-domain  into  the  x -domain.  These  Fourier  conjugate  domains  lie  in  the  axial 
direction  of  the  cylinder.  The  phenomena  of  special  concern  in  the  present  report  are 
those  relating  to  aliasing  and  pass  and  stop  bands  in  that  modal  response  [2-6].  Hybrid 
and  natural  cylinders  are  exemphfied.  A  hybrid  cylinder  is  defined  by  analytically 
removing  the  membrane  free  waves  from  the  cylinder.  In  this  sense,  the  mechanical 
response  of  a  hybrid  cylinder  is  akin  to  that  of  a  panel  that  is  enjoying  the  same 
plating  [2].  This  report  assumes  that  the  descriptions  of  these  phenomena  and  the 
analytical  models  of  the  cylindrical  structures  are  familiar  to  the  reader;  nonetheless  a 
cursory  analytical  development  is  included.  A  more  detailed  analytical  development  can 
be  found  in  References  2-6.  However,  the  analytical  development  of  the  double-sum 
technique  that  is  presented  in  Reference  1  is  not  repeated  here;  Reference  1  is 
considered  an  essential  companion  to  this  report  and  the  reference  to  equations  in 
Reference  1  are  prefaced  by  unity.  Thus,  Eq.  (6)  in  Reference  1  is  designated  Eq.  (1-6)  in 
the  present  report.  The  physical  construction  of  the  externally  driven  ribbed  cylinder 
and  the  chosen  coordinate  system  are  sketched  and  depicted  in  Fig.  1,  respectively. 

In  Section  II  the  response  and  the  external  drive  that  generate  it  are  expanded  in 
terms  of  circumferential  eigenfunctions  [7].  The  coefficients  in  these  expansions  are  the 
modal  responses  and  modal  external  drives.  The  modal  properties  of  the  cylindrical  shell 
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that  relate  the  modal  response  to  the  corresponding  modal  external  drive  are  defined; 
this  relationship  defines  the  modal  impulse  response  function.  It  is  argued  that  for  a 
regularly  ribbed  cylinder,  the  expressions  for  the  impulse  response  function  and  its  terms 
and  factors  are  readily  stated  and  interpreted  in  the  [k,  n,  co}  -domain,  where  (to)  is  the 
frequency  variable,  (n)  is  the  mode  index  and  (k)  is  the  axial  wavenumber  variable  [8]. 
Indeed,  in  this  domain  the  relationship  between  the  modal  response  (k,  co)  and  the 
modal  external  drive  (k,  co)  may  be  expressed  in  the  form 

V„(k,co)  =  G„„(k,co)  [P^(k,co)  -  P,,ik,co)]  , 

where  P^{k,co)  is  the  drive  in  lieu  of  the  ribs  and  G^^(k,co)  is  the  modal  admittance  of 
the  uniform  (unribbed)  cylinder.  The  external  drive  is  adequately  defined  in  the  simple 
form 


Pen  (k,C0)  =  P^n  ( ® )  exp{ikx^ )  , 

where  (x^ )  is  the  position,  in  the  x  —  domain,  to  which  the  external  drive  is  applied  and 
it  is  assumed  that  the  modal  strength  Pen((0)  is  independent  of  (k).  The  uniform 
cylinder  modal  admittance  G^n(^^^)  comprises  two  distinct  terms;  the  modal 
mechanical  surface  impedance  Zpn{k,co)  and  the  modal  fluid  surface  impedance 
Zl  (k,  co)  on  the  outside  surface  of  the  cylinder 

G„(k,0))  =  [Zp„ (*,<»)+  Zf(k,<t))]-'  . 
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On  the  inside  of  the  cylinder  vacuum  is  assumed  to  prevail.  The  modal  mechanical 
surface  impedance  Zp^(k, co),  in  turn,  also  comprises  two  terms 

Zp,{k,co)  =  ZP(k,(0)  +  z;:(k,Q}) , 

where  Z^{k,(o)  is  the  modal  flexural  surface  impedance  and  Z^{k,(o)  is  the 
membrane  surface  impedance  [8-10].  When  the  membrane  surface  impedance  is 
analytically  removed;  (Jc,  (o)  =  0 ,  the  admittance  (k,  co),  given  now  by 

G„(k,o>)  =  +  Zf(k,0»r'  , 

describes  the  fluid  loaded  modal  admittance  of  a  uniform  (unribbed)  hybrid  cylinder. 

The  drive  co)  in  lieu  of  the  ribs  is  not  only  a  functional  of  the  properties  of 
the  uniform  cylinder;  defined  by  G„^(k,co),  but  is  naturally  also  a  functional  of  the 
properties  of  the  ribs.  In  a  regularly  ribbed  cylinder  the  properties  of  the  ribs  are  defined 
in  terms  of  the  modal  line  impedance  (co)  and  the  separation  distance  (b)  between 
adjacent  ribs.  This  line  impedance  can  be  either  mass,  resistance  or  stiffness  controlled; 
the  change  from  one  characteristic  to  another  is  frequency  dependent  [11].  The  modal 
drive  P^(k,(o),  in  addition,  is  a  functional  of  the  external  drive.  In  the  definition  of  the 
external  drive  here  proposed,  the  dependence  on  the  position  of  application  (x^ )  is  of 
particular  significance. 

It  is  convenient  to  normalize  the  modal  response  in  the  form 

V,(k,co)  =  V^,,(k,(o)  -  V,,(k,(o)  =  G,„(/:,£o)  {P,,(k,co)  -  P,,(k,co)]  ; 

V^n(k,co)  =  G^^{k,(0)  Penile, (0)  ]  V,„(k,(o)  =  G„^ik,(0)  Penile, co)  . 
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This  normalization  of  the  modal  response  is  achieved  by  the  appropriate  normalization  of 
its  constituent  terms  and  factors,  namely 

V,{K(o)  =  [V,{K(o)IV,,{(0)\  ;  V^io))  =  [PMKicom)]  ; 

G^(k,<o)  =  [(i(om)  G^,{k,(o)]  ;  Uk,(o)  =  [P^ik,co)/ PeW]  , 

where  (to  m)  is  the  surface  mass  impedance  of  the  plating  of  the  cylinder  and  Pf„(cD),  as 
already  mentioned,  is  the  strength  of  the  modal  external  drive. 

The  absolute  values  of  V„„(k,Q)),  or  equivalently  G^(k,Q)),  of  P^^{k,(0)  and  of 
V^(k,co)  are  computed  and  displayed  as  functions  of  (ak)  in  a  frequency  waterfall 
format,  where  (a)  is  the  radius  of  the  cylinder  [12].  The  frequency  (co)  is  normalized  by 
the  critical  frequency  (o^ )  of  the  flexural  free  waves  with  respect  to  the  speed  of  sound 
(c)  in  the  fluid.  The  discussions  of  these  and  other  presentations  in  this  report  are 
limited  to  the  frequency  range  defined  by 

{2(0 J(0^)  <{(0l(0^)  <  0.6  ;  (o,  =  {cjd)  , 

and  to  the  wavenumber  range  defined  by 

0  <:{ak)  <,  75  , 

where  co^  is  the  ring  frequency  of  the  cylinder.  The  speed  q  =  [Yh! , 
where  Y ,v,  m  and  h  are  the  Young’s  modulus,  the  Poisson’s  ratio,  the  surface  mass  and 
the  thickness  of  the  plating  of  the  cylinder,  respectively,  is  the  longitudinal  speed  in  that 
plating.  One  may  find  it  convenient  to  list  a  set  of  standard  parametric  values  to  define 
standard  displays.  When  variations  in  the  parametric  values  are  instituted,  it  is  necessary 
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to  report  only  those  parameters  that  are  varied.  The  remaining  parameters  are  assumed 
to  maintain  standard  values. 

The  normalized  modal  response  {k,  (o)  of  a  uniform  (unribbed)  cylinder  under 
standard  parametric  values,  and  a  few  specific  variations  on  these  parametric  values,  are 
presented  in  Figs.  2  and  3.  The  first  set;  i.e..  Fig.  2,  pertains  to  a  hybrid  cylinder  and  the 
second  set;  i.e..  Fig.  3,  pertains  to  a  natural  cylinder.  The  footprints  of  fluid  loading  and 
membrane  free  waves  are  of  particular  interest  in  these  two  sets  of  figures. 

The  normalized  modal  drive  {k,  (o)  in  lieu  of  the  ribs  under  standard  parametric 
values,  and  a  few  specific  variations  on  these  parametric  values,  are  presented  in  Figs.  4 
and  5.  The  third  set;  i.e..  Fig.  4,  pertains  to  a  hybrid  cylinder  and  the  fourth  set;  i.e.. 
Fig.  5,  pertains  to  a  natural  cylinder.  The  footprints  of  aliasing  and  stop  and  pass  bands 
are  of  particular  interest  in  these  two  figures. 

The  normalized  modal  response  V^{k,(0)  of  the  regularly  ribbed  cylinder  under 
standard  parametric  values,  and  a  few  specific  variations  on  these  parametric  values,  are 
presented  in  Figs.  6  and  7.  The  fifth  set;  i.e..  Fig.  6,  pertains  to  a  hybrid  cylinder  and  the 
sixth  set;  i.e..  Fig.  7,  pertains  to  a  natural  cylinder.  The  footprints  of  the  normalized  drive 
Pg^(k,co),  as  depicted  in  Figs.  4  and  5,  are  clearly  discernible  in  the  normalized  response 
V„(k,co).  However,  the  strict  aliasing  in  the  drive  is  impaired  in  the  response.  This 
defect  is  perpetrated  by  the  term  that  describes  the  normalized  modal  response 
Kon(^>®)  ill  til®  absence  of  ribs  and  by  the  factor  G„^(k,(o)  that  multiplies  the 
normalized  aliased  modal  drive  P^^(k,Q))',  neither  the  term  nor  the  factor  are  aliased 
quantities.  Does  the  impairment  of  the  aliasing  in  (k,  co)  imply  a  corresponding  defect 
in  the  pass  and  stop  bands  of  the  modal  response  as  compared  with  that  exhibited  by 
the  modal  drive  in  lieu  of  the  ribs?  The  answer  to  this  question  is  sought,  in  part,  by 
using  complementarity  to  Figs.  2  through  7.  This  complementarity  is  derived  by  a 
Fourier  transformation  of  V^^(k,0})  =  [G^„(k,Q})  P^^(k,co)],  P^^{k,(o)  and  V„{k,co) 
into  v^„(x,co),  Pg„(x,co)  and  v^(x,(o),  respectively.  The  Fourier  transformation  of 
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these  quantities  is  carried  out  in  terms  of  the  double-sum  technique  developed  in 
Reference  1. 

In  Section  III  the  procedure  for  computing  by  a  Fourier  transformation 

of  V„^(kyCo),is  explained  and  the  results  are  displayed  in  Figs.  8  through  10.  Several 
variations  in  the  parametric  values  are  investigated  and  comparisons  between  Figs.  2 
and  3  that  feature  V^,^(k,co)  and  the  corresponding  Figs.  8  through  10  that  feature 
v^„(jc,£o)  are  conducted  and  discussed.  Interferences  between  flexural,  fluid  and 
membrane  free  waves  at  regions  where  their  amplitudes  are  of  comparable  strengths  are 
of  particular  interest.  Such  interferences  are  revealed  when  the  dominant  flexural  free 
waves  are  subdued  by  selectively  increasing  the  damping  that  is  assigned  to  them.  It  is 
thus  found  that  a  complemental  investigation  of  V^^{k,(o)  and  v^^{x,(o)  is  an 
analytical  tool  of  value. 

In  Section  IV  the  procedure  for  computing  (x,  (o)  by  a  Fourier  transformation 
of  (O),  is  explained  and  the  results  are  displayed  in  Figs.  4,  5,  11  and  12.  In  Figs.  4 
and  5  computations  of  p^^(x,co)  are  superimposed,  frequency- wise,  on  the 
corresponding  figures  that  display  P^^{k,co).  Figures  4  and  5  exhibit  clearly  the 
phenomena  of  aliasing  and  pass  and  stop  bands.  The  complementarity  between  the 
Fourier  conjugate  pairs,  P^^(Ji,(o),  on  the  one  hand,  and  p^^(x,(o),  on  the  other,  are 
more  directly  and  simultaneously  contrasted  in  these  figures.  The  value  of 
complementarity  as  an  analytical  tool  is  thus  supported  in  Figs.  4  and  5.  The  influence 
of  fluid  loading  on  the  phenomenon  of  pass  and  stop  bands  is  investigated  in  Figs.  1 1 
and  12.  Moreover,  the  bottoming  of  the  flexural  stop  bands  by  the  fluid  and  membrane 
free  waves  are  discernible  in  these  figures. 

In  Section  V  the  procedure  for  computing  v„(x,  (o)  by  a  Fourier  transformation  of 
V,^(k,  co),  is  explained  and  the  results  are  displayed  in  Figs.  13  through  15.  Of  particular 
interest  is  the  bottoming  of  the  response  in  a  flexural  stop  band  by  the  fluid  and 
membrane  free  waves.  The  similarity  of  this  bottoming  to  the  corresponding  bottoming 
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induced  by  damping  the  flexural  free  waves,  is  striking,  [cf.  Figs.  8  through  10  and  13 
through  15.]  The  complementarity  between  Figs.  6  and  7  and  the  corresponding 
Figs.  13  through  15,  again  shows  the  value  of  its  employment.  Of  particular  interest  in 
this  final  section  are  the  comparisons  between  Figs.  4  and  5  and  11  and  12,  on  the  one 
hand,  and  Figs.  6  and  7  and  13  through  15,  on  the  other.  These  comparisons  reveal  that 
the  impairment  of  the  aliasing  patterns  in  (k,  co) ,  as  compared  with  those  in  (k,  co) , 
is  not  duplicated  in  the  pass  and  stop  bands  in  v^(x,co),  as  compared  with  those  in 
Pj.^(XyCO).  Indeed,  if  at  all,  the  phenomenon  of  pass  and  stop  bands  in  v„(x,a))  is 
sharper  than  it  is  in  (jc,  co) . 
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n.  MODAL  RESPONSE  OF  A  RIBBED  CYLINDER 
The  circumferential  uniformity  and  finiteness  of  the  cylinder,  unribbed  and  ribbed, 
allow  quantities  to  be  expanded  in  terms  of  circumferential  eigenfunctions  where 

0  is  the  circumferential  variable  as  shown  in  Fig.  1.  Thus,  the  response  v(x,(f),Q))  and  the 
external  drive  Pg{x,^,(0)  can  be  modalized  in  the  x- domain 

v(x,^,£0)  =  V„(x,6))  ; 

n 

p^ix,(l),0))  =  2  Pen(x,CO)  ,  (la) 

n 

and  in  the  k  -domain 

V^k,^,o>)  =  2  <!>«(«  v„(k,co)  ; 


PAk,(tl,<0)  =  2  Pen(k.(»)  . 

n 

Quantities  in  Eq,  (1)  are  related  Fourier  transforms;  in  particular 

(lb) 

V„(X,C!))  =  (271)  j 

dk  exp(-ixk)  V„(k,(jo)  , 

(2a) 

Pen(x,(a)  =  (27t)~^'‘^ 

J  dk  exp(-ixk)  Pg^(k,co)  . 

(2b) 
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The  circumferential  eigenfunctions  O„(0),  are  assumed  to  be  orthogonal  and  complete 


I  <!>„(«  #  <pU<I»  = 


0  <  (^-f)  <  2;r  , 


(3) 


where  5^  and  5  are  the  Kronecker  delta  and  the  delta  function,  respectively  [8].  The 
modal  response  v^{x,(o)  of  the  cylinder  to  the  modal  external  drive  Pg^(x,co)  in  the 
X  -  domain  can  be  formally  expressed  in  the  form 


g^(x\x',CO)  dx'  Pen(x\CO)  , 


(4a) 


and  in  the  k  -domain  in  the  form 


v„(k.0)  =  jG„(k\k',0))  dk'  p„(k',0))  .  (4b) 

where,  in  addition  to  Eq.  (2),  the  impulse  response  functions  in  these  equations  are 
related  Fourier  transforms;  e.g,, 

=  (27r)"^  J  J  ^xp[-i(xk - x'k')]  Gn(k\k',(o)  .  (5) 
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It  is  convenient  to  normalize  the  modal  impulse  response  function  by  the  surface  mass 
impedance  {icom)  of  the  plating  of  the  cylinder  and  the  modal  drive  by  the  modal  ring 
drive  namely 

G^{,k\k\(o)  =  {i(om)  G^{,k\k\(o)  ; 

P,,{k,co)  =  [P,„a,a))/F;„(ft>)]  =  exp{ix,k)  ,  (6a) 

where  (w)  is  the  surface  mass  of  the  plating  and  )  is  the  position  of  the  ring  to  which 
the  modal  external  drive  is  applied.  The  normalization  of  the  modal  response  V^(k,co) 
follows 

V,(k,co)  =  [V,(k,co)/V^(co)]  ;  V^(co)  =  [Pf,(co)/(io}m)]  .  (6b) 

The  modal  impulse  response  function  G^(k\  k',0))  for  a  regularly  ribbed  cylinder  is 
readily  derived  in  the  k  -domain  in  the  simple  operator  form 

G,(k\k',0))  =  G^,{k,(o)  [1-5, (^)  {/,„(^,(»)  (/:,«)]  5{k-k')  ,  (7) 

where 


=  [z^mib]  G.„(k,(o) , 

(8a) 

hnik,(0) 

=  [l+Si,(k)  [Q.„(k,(o}]]-'  , 

(8b) 

the  wavenumber  operator  5,(^)  is  defined  in  Eq.  (1-3),  {b)  is  the  separation  distance 
between  adjacent  ribs,  Z^„(£t})  is  the  modal  ring  impedance  of  the  ribs,  and  G^^{k,  co)  is 
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the  normalized  modal  surface  admittance  of  the  uniform  (unribbed)  cylinder.  The 
incomplete  curly  brackets  preceded  by  the  wavenumber  operator  Sj,  (k)  in  Eq.  (7) 
merely  indicates  that  the  wavenumber  operation  to  be  performed  is  yet  to  be  completed; 
indeed,  G„(^  I  ^',0))  is  a  wavenumber  operator  and  not  just  an  algebraic  factor.  The 
normalized  quantity  (o)  is  the  inverse  of  three  terms 

G^^{K(o)  =  {i(om)  G^„(k,co)  = 

[Z^(k,a)  +  Z“(i,<B)  +  Z/ (A:, <»)]-'  ,  (9a) 

where  Z^(k,a>),  Z^(k,Q})  and  Z/  (k,oj)  are  the  normalized  flexural,  membrane  and 
fluid  modal  surface  impedances,  respectively,  of  the  uniform  cylinder.  A  hybrid  cylinder 
is  defined  by  a  (natural)  cylinder  from  which  the  membrane  free  waves  are  analytically 
removed  by  setting  (k,  co)  equal  to  zero  in  Eq.  (9a).  A  hybrid  cylinder  is  defined  by 
a  normalized  modal  surface  admittance 

=  iz;(k,<0)  +  2/(1:, n))]-'  ;  Z/d:,®)  =  0  .  (9h) 

Recognizing  that  the  mechanical  surface  impedance  Zp(k,ky,(o)  for  a  uniform  panel 
relates  to  (k,  co)  in  the  form 

Zp(k,(nla),co)  ^  ZP{k,co)  ,  (10) 

and  that  the  fluid  surface  impedance  Zj-{k,(n  I  a),co)  on  a  panel  qualitatively  relates  to 
z/  (fc,6)),it  is  inferred  that  a  hybrid  cylinder  is  akin  to  a  panel,  and  in  some  cases,  the 
results  obtained  on  a  hybrid  cylinder  can  be  inferentially  translated  to  those  on  a  panel 
[2].  The  quantity  Q^n(^,co)  is  the  ratio  of  the  equivalent  modal  surface  impedance 
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[Z^{(o)lb\  of  a  rib  to  the  modal  surface  impedance  ^  of  the  uniform 

cylinder.  It  is  convenient  to  express  in  the  form 


'Z‘sni(o)  =  (icoM)  A^icolcoJ  ,  (11) 

and  hence  from  Eq.  (8a)  one  derives 

{Ml mb)  A^icofaJ  ,  (8c) 

where  M  is  a  line  mass  and  A^{(ol (O^)  is  a  resonance  bearing  factor.  In  this  report 
A„(co/(»„)  is  equivalently  set  either  to  unity  or  to  -[l  +  (a„ct)/(»„)^]“\  the  ring 
impedance  of  the  ribs  in  the  former  setting  is  mass  controlled  and  in  the  latter  stiffness 
controlled  [2].  Substituting  Eq.  (7)  in  Eq.  (4b)  one  derives 

V^ik,(o)  =  V„^{k,(o)  -  V^^{k,co)  ,  (12) 

V^^(k,co)  =  G^„(k,(o)  exp(ix^k)  ,  (13a) 

V,^(k,(o)  =  G^^{k,(o)  P^{k,(o)  ,  (13b) 

where 

p,,{k,(o)  =  [p,,{k,(o)ip:nm  = 

!2oon(^,ty)  expiix^k)]  ,  (14a) 

is  the  normalized  modal  drive  in  lieu  of  the  ribs.  Clearly,  P^^{k,(o)  is  an  aliased  quantity 
in  {k)  with  respect  to  (k^)-,  =  (In lb),  [cf.  Eq.  (l-17a).]  Moreover,  one  recognizes 
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that  If,^ik,(o)  as  stated  in  Eq.  (8b)  is  similarly  aliased  in  {k)  with  respect  (kTi). 
[cf.  Eq.  (l-19a).]  In  consequence  P^(k,a>)  can  be  factorized  in  the  form 

=  hn^k,(0)  Si,{k)  [Q^rtik,(0)  exp{ix^k)]  .  (14b) 

[cf.  Eq.  (l-18a).]  It  is  apparent  from  Eqs.  (12)  through  (14)  that  the  behavior  of  the 
modal  response  V^{k,(o)  and  a  number  of  terms  and  factors  in  its  composition;  e.g.,  the 
modal  response  V^^(kyCO)  in  the  absence  of  ribs  and  the  drive  Pj^(k,(0)  in  lieu  of  the 
ribs,  are  merely  functional  of  the  properties  of  the  fluid  loaded  uniform  cylinder,  the  ribs 
and  the  external  drive.  These  properties  are  described  by  the  modal  surface  admittance 
G^J<k,(o)  of  the  uniform  cylinder,  the  modal  line  impedance  Z^„(to)  and  the  separation 
distance  {b)  of  the  ribs  and  the  ring  position  {x^ )  of  the  external  drive.  The  absolute 
values  of  the  quantities  V^^{k,(o)  or  equivalently  G„^(k,co),  P^^(k,(o)  and  V^(k,(o),SLS 
functions  of  (ak),  are  computed  and  presented  in  typical  frequency  waterfall  displays  in 
Figs.  2  and  3,  4  and  5,  and  6  and  7,  respectively  [12].  The  even  numbered  figures  in  each 
set;  e.g..  Fig.  2,  relate  to  a  hybrid  cylinder  and  the  odd  numbered;  e.g..  Fig.  3,  relate  to  a 
natural  cylinder.  A  typical  case  is  defined  by  standard  parametric  values.  These 
standard  values  may  also  serve  to  define  other  cases  by  merely  selecting  changes  in 
specific  standard  values.  The  standard  parametric  values  and  their  definitions 
follow:  the  flexural  free  waves  are  defined  by  a  wavenumber  kp, 

=  kpo(^-ilp) --{(ocojc^); 

{hid)  =  10"^=  rip  =  10"^  ,  (15a) 
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the  longitudinal  free  waves  are  defined  by  a  wavenumber 


k,  =  k^(l-iVt)  ;  k,^  =  (ffl/c,)  :  c?  =  (r/p.Xl-v^)-'  ; 

m  =  /iPj  ;  (c^lc)  =  3.5  ;  =  5  x  10^  ,  (15b) 

the  shear  free  waves  are  defined  by  a  wavenumber 

=  kso(^-Ms)>  Ko  =(0}lc,);  (cjc^f  =  [(l-v)/2]  ; 

V  =  (1/3)  ;  7?,  =  5  X  10"^  ,  (15c) 

the  fluid  loading  is  defined  by  a  fluid  loading  parameter 

=  (pclco^m)  ;  {a(oJc)  =  97.53  ;  =  10"^  ,  (15d) 

the  ring  frequency  (Oj.  and  the  critical  frequency  (O^  are  defined 
(O^  ={cja)\  (o^  =  Vl2  (c^ //ic^)  ; 

{(0j(0^)  =  3.6  X  10'^  ,  (15e) 

and  the  definition  of  the  ribs  and  mode  index  {n)  are  specified 

{bid)  =  0.3  ;  {Ml mb)  =  0.2  ;  {xjb)  =  0.5  ; 

A^{(0l(0^)  =  1  ;  n  =  1  ,  (15f) 
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where  (a)  is  the  radius,  Y  is  Young’s  modulus,  (v)  is  the  Poisson’s  ratio  and  (p^)  is  the 
density  of  the  plating  of  the  cylinder  and  (p)  and  (c)  are  the  density  and  sound  speed  in 
the  fluid.  The  interest  in  this  report  is  confined  to  the  spectral  range 

{0  ,  (2(0, 1(0,)]  <  {(ak)  ,  ((ol(o,)}  <  {75  ,  0.6}  ,  (15g) 

notwithstanding  that  some  displays  are  extended,  but  not  discussed,  into  the  frequency 
range  ((0l(0,)  <(2(0,l(0,).  Computations  and  their  representations  in  this  low 
frequency  are  discussed  under  a  separate  cover  [2].  The  first  sub-figures  in  the  first  and 
second  sets;  i.e..  Figs.  2a  and  3a,  are  computed  under  standard  parametric  values.  The 
flexural  free  waves  and  the  fluid  free  waves  are  clearly  discernible  in  Fig.  2a.  In  Fig.  3a, 
in  addition  to  the  flexural  and  fluid  free  waves,  membrane  free  waves  are  clearly 
discernible.  The  membrane  fi'ee  waves  comprise  the  longitudinal  and  shear  free  waves. 
“Curvature  free  waves”  that  supersede  the  flexural  free  waves  and  the  origins  and  the 
higher  surface  impedances  of  the  membrane  free  waves  are  found  largely  to  reside,  in 
Fig.  3a,  in  the  low  frequency  range;  (col  (O,)  <  (2(o,  I  (O,).  Again,  discussion  of  these 
features  are  not  conducted  in  this  report  [2].  Figures  2a  and  3a  are  repeated  in  Figs.  2b 
and  c  and  3b  and  c,  respectively,  except  that  the  fluid  loading  parameter  (e, )  is  changed 
from  the  standard  value  of  10”^  to  10”^  and  10“\  respectively.  In  this  connection  it  is 
observed  that  a  fluid  loading  with  e,  =  10“^  is  light,  with  e,  =  10~^  is  moderate  and 
with  e,  =  10"^  is  heavy.  The  change  of  from  10'^  to  lO”^  renders  the  fluid  free 
waves  hardly  discernible  in  Figs.  2b  and  3b,  little  else  is  different  between  these  figures 
and  Figs.  2a  and  3a,  respectively.  The  change  of  (e^)  from  10“^  to  10“^  intensifies  the 
fluid  free  waves  in  Figs.  2c  and  3c.  Moreover,  this  change  in  fluid  loading  shifts  the  free 
waves,  especially  the  subsonic  flexural  waves,  to  a  higher  wavenumber  range  and 
weakens  the  strengths  of  the  supersonic  membrane  free  waves.  The  shifts  are  due  to  the 
“added  surface  mass”  and  the  weakening  is  due  to  the  “radiation  damping”  effects 
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associated  with  the  fluid  loading.  Figures  2d  and  3d  repeat  Figs.  2a  and  3a,  respectively, 

_ <1 

except  that  the  flexural  loss  factor  (rip )  is  changed  from  the  standard  value  of  10  to 
3  X  10“^;  a  substantial  increase  in  the  damping  of  the  flexural  free  waves.  Indeed,  the 
flexural  free  waves  in  Figs.  2d  and  3d  are  substantially  subdued  compared  with  those  in 
Figs.  2a  and  3a.  Figures  2a  and  3a  are  repeated  in  Figs.  2e  and  3e,  respectively,  except 
that  the  mode  index  (n)  is  changed  from  the  standard  value  of  unity  to  zero.  This 
change  hardly  induces  any  change  in  the  flexural  free  waves;  however,  a  drastic  change 
is  induced  in  the  membrane  free  waves.  The  membrane  free  waves  in  Fig.  3a  are 
comprised  of  longitudinal  and  shear  free  waves,  in  Fig.  3e  of  longitudinal  free  waves 
only.  The  shear  free  waves  join  the  longitudinal  when  the  mode  index  (n)  equals  to  or 
exceeds  unity.  Thus,  in  Fig.  3f,  which  repeats  Fig.  3a  except  that  the  mode  index  (n)  is 
set  equal  to  (2),  the  membrane  free  waves  are  composed  of  both,  the  longitudinal  and 
shear  free  waves.  Finally,  the  first  set  is  completed  by  Fig.  3g.  This  figure  repeats 
Fig.  3e,  except  that  the  flexural  loss  factor  (77^  )  is  changed  from  the  standard  value  of 
10  to  3  X  10  .  Again,  the  flexural  free  waves  are  substantially  subdued  by  this 
change,  [cf.  Figs.  2d  and  3d.] 

The  first  sub-figures  in  the  third  and  fourth  sets;  i.e..  Figs.  4a  and  5a,  are  computed 
under  standard  parametric  values.  The  aliasing  of  the  normalized  modal  drive  P^^(k,co) 
in  lieu  of  the  ribs  is  apparent  in  these  figures;  the  dominant  patterns  are  those  associated 
with  the  flexural  free  waves.  Nonetheless,  aliasing  in  the  membrane  free  waves  are 
discernible  in  Fig.  5a.  Aliasing  is  basically  a  scattering  phenomenon  and  although  the 
aliasing  patterns  of  the  membrane  free  waves  are  scant,  the  scattering  of  the  membrane 
free  waves  by  “simple”  ribs  are  clearly  established  in  Fig.  5  a.  In  part,  the  visibility  of  the 
membrane  free  waves  in  Fig.  5a  is  aided  by  the  locus  minima  in  P^^(k,co).  These  locus 
minima  occur  in  Pg^(k,co),  in  both  Figs.  4a  and  5a,  whenever 
(ak)  =  ±  (aKi)  [(2j+l)/2]  and  is  due  to  the  symmetry  in  this  quantity.  The 
symmetry  results  from  the  standard  position  of  the  external  drive  at  (x^  lb)  =  0.5. 


17 


This  symmetry  is  removed  in  Figs.  4b  and  5b.  These  figures  repeat  Figs.  4a  and  5a, 
respectively,  except  that  (x^lb)  is  changed  from  0.5  to  0.3.  The  locus  of  minima  are 
absent  in  the  patterns  of  Figs.  4b  and  5b.  Moreover,  the  patterns  relating  to  the 
membrane  free  waves  in  Fig.  5b  are  weakened.  The  fluid  free  waves  are  hardly 
discernible  in  Figs.  4a  and  b  and  5a  and  b.  A  major  feature  in  the  patterns  of  these 
figures  is  the  cancellations  between  the  peaks  of  a  pair  of  flexural  aliasing  orders  in  the 
region  in  which  they  cross.  Regions  of  these  kind  lie  regularly  in  certain  frequency 
bands.  The  cancellations  are  made  possible  by  the  interactions  of  the  ribs  among 
themselves  via  the  cylindrical  shells;  in  the  absence  of  such  interactions  the  cancellations 
are  also  absent  [2,13].  Figures  4a  and  5a  are  repeated  in  Figs.  4c  and  5c  respectively, 
except  that  the  fluid  loading  parameter  is  increased  from  the  standard  value  of  10“^  to 
10“^  The  influence  of  this  increase  is  largely  discernible  in  the  shifts  in  the  frequency 
bands  at  which  pairs  of  various  flexural  aliasing  orders  cross;  the  shifts  are  to  lower 
frequencies  and  are  consistent  with  the  wavenumber  shifts  observed  between  Figs.  2a 
and  3a  and  Figs.  2c  and  3c.  In  addition,  the  flexural  free  waves  in  both  figures  and  the 
membrane  free  waves  in  Fig.  5c  are  subdued  by  this  increase  in  the  fluid  loading. 
Remembering  that  under  the  influence  of  ribs  the  subsonic  flexural  free  waves  can 
become  reasonable  radiators,  the  subdued  flexural  and,  again,  the  membrane  free  waves 
can  be  attributed  to  the  increase  in  radiation  damping  with  the  increase  in  fluid  loading. 
However,  the  fluid  free  waves,  as  such,  are  stiU  hardly  discernible  in  Figs.  4c  and  5c. 
[cf.  Figs.  4a  and  b  and  5a  and  b,  respectively.]  Figures  4d  and  5d  repeat  Figs.  4a  and  5a, 
except  that  the  standard  mass  controlled  ring  impedance  of  the  ribs;  i.e.,  with 
A^(colco^)  =  1,  is  changed  to  that  of  stiffness  controlled;  i.e.,  with 
A^{(ol(0^)  -  -  [1  +  (5fi)  /  n)p )  ]“  .  The  change  is  largely  discernible  in  the  shifts  in  the 
frequency  bands  at  which  pairs  of  various  flexural  aliasing  orders  cross;  the  shifts  are  to 
higher  frequencies.  Otherwise  the  influence  of  the  change  on  patterns  of  the  aliased 
modal  drive  P^{k,(o)  is  insignificant,  notwithstanding  that  the  cancellations  are  weaker 
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and  narrower  at  the  crossings  of  pairs  of  various  flexural  aliasing  orders  as  the  ring 
impedance  of  the  ribs  is  diminished,  by  definition,  at  the  higher  frequency  ranges. 
Figure  5a  is  repeated  in  Fig.  5e,  except  that  the  mode  index  («)  is  changed  from  the 
standard  value  of  unity  to  zero.  The  absence  of  the  shear  free  waves  in  Fig.  5e  is 
consistent  with  Fig.  3e. 

The  first  sub-figures  in  the  fifth  and  sixth  sets;  i.e..  Figs.  6a  and  7a,  are  computed 
under  standard  parametric  values.  The  patterns  in  the  normalized  modal  response 
V„(k,(o)  manifest  footprints  of  the  aliasing  in  the  normalized  drive  P^„(k,(o).  However, 
since  neither  the  term  P^^(k,(o)  [=exp{ix^k)]  nor  the  factor  G^„(A:,a}),in  Eqs.  (13)  and 
(14),  are  aliased  in  (k)  with  respect  to  (k^),  the  aliasing  in  P^(k,0})  is  spoiled  in 
V^(k,(o)’,  the  quantity  V^{k,(o)  is  not  strictly  aliased.  The  “ribs’  free  term”  V^^ik,0))  in 
this  quantity  accentuates  and  reemphasizes  the  zeroth  order  patterns  and  the  factor 
G^„(k,co)  in  the  “ribs’  term”  V^^{k,co)  subdues  the  aliased  patterns  in  P^{k,(o)  once 
[{ak)^  -t-  (n/a)^]  >  (kp^)^  is  satisfied.  The  subduing  increases  the  more  this 
inequality  is  satisfied.  The  first  of  these  effects  is  manifested  in  that  the  patterns  in 
Figs.  2a  and  3a,  are,  respectively,  overlaid  in  Figs.  4a  and  5a.  The  second  is  manifested  in 
that  the  patterns  in  Figs.  4a  and  5a  are  suppressed  in  Figs.  6a  and  7a  in  spectral  ranges 
that  lie  beyond  the  flexural  free  waves.  The  second  effect  is  made  clearer  in  Figs.  6a.  1 
and  7a.  1.  These  figures  are  the  clipped  version  of  Figs.  6a  and  7a,  respectively;  values 
that  lie  below  a  given  threshold  in  Figs.  6a  and  7a  are  removed  in  Figs.  6a.  1  and  7a.  1 
[2,3].  For  the  sake  of  completeness  Figs.  6a  and  7a  are  repeated  in  Figs.  6b-d  and  7b-e, 
with  changes  that  correspond  to  those  made  in  Figs.  4b-d  and  5b-e  with  respect  to 
Figs.  4a  and  5a,  respectively.  The  features  in  the  patterns  of  Figs.  6a  and  7a  drawn  from 
the  patterns  in  Figs.  2a  and  4a  and  3a  and  5a,  respectively,  can  be  generalized  to  the 
other  sub-figures  in  these  categories;  e.g.,  the  patterns  in  Fig.  7c  are  overlaid  of  patterns 
in  Fig.  3c  and  of  patterns  in  Fig.  5c.  The  patterns  in  the  latter  figure  are,  however, 
subdued  in  spectral  ranges  that  lie  beyond  the  flexural  free  waves.  Thus,  changes  in  the 


19 


parametric  values  and  accompanied  influences  in  Figs.  2  through  5,  if  correspondingly 
selected,  impose  similar  influences  in  Figs.  6  and  7. 

The  derivation  and  manipulation  of  the  modal  response  V^{k,(0)  and  its 
composition  are  found  to  be  relatively  simple  in  the  k  -domain.  The  algebraic  nature  of 
the  modal  response  V^^{k,co)  induced  by  the  normalized  modal  external  drive  Pg^{k,0)) 
in  the  absence  of  the  ribs,  the  harmonic  nature  of  the  drive  P^^{k,(o)  in  lieu  of  the  ribs 
and,  again,  the  algebraic  nature  of  the  modal  response  V^^(k,co)  induced  by  this  drive, 
render  efficacious  the  ^-domain  for  formulating  the  modal  response  V^{k,0}).  The 
descriptions  of  this  quantity  and  terms  and  factors  in  its  composition  are  presented  in 
Figs.  2  through  7.  These  descriptions  demonstrate  that  the  flexural  free  waves  are  the 
dominant  response  contributors,  that  the  fluid  free  waves  are  weak  contributors  even 
when  fluid  loading  is  heavy  and,  finally,  although  the  membrane  free  waves  are 
discernible  and  they  are  evidently  scattered  by  the  ribs,  the  membrane  free  waves  are 
nonetheless  weak  contributors  to  the  response,  especially  at  the  higher  frequency  range, 
where  (2co^  /®c)  5(®/®c)-  Figs.  2f  and  3f  and  g  the  flexural  free  waves  are 
substantially  weakened  by  increasing  the  flexural  loss  factor  (rip)  by  more  than  an 
order  of  magnitude.  The  other  free  waves  remain  substantially  unaltered  so  that  they  are 
relatively  accentuated.  The  quantitative  assessments  of  the  weakening  of  the  flexural 
free  waves  and  the  accentuation,  thereby,  of  the  other  free  waves  are  not  readily 
afforded  by  Figs.  2  and  3  which  include  Figs.  2f  and  3f  and  g.  In  this  light,  can  a  Fourier 
transformation  of  into  the  x- domain  assist  in  these  assessments? 
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m.  MODAL  RESPONSE  OF  A  UNIFORM  CYLINDER  IN  THE  SPATIAL  DOMAIN 

The  expression  for  the  normalized  modal  response  V„^{k,(o)  of  a  uniform  cylinder 
is  stated  in  Eq.  (13a).  As  already  observed,  neither  the  term  exp(ix^k)  nor  the  factor 
G„^{k,(o)  are  aliased  and,  therefore,  the  modal  response  V^(k,(o)  in  the  -domain  is 
not  aliased.  The  Fourier  transformation  of  V^{k,(o)  into  the  spatial  x -domain  affords 
no  factorial  saving  performing  the  transformation  either  via  Eq.  (1-lOb)  or  via  Eq.  (1-14); 
indeed,  computations  of  this  kind  are  relatively  simple.  Nonetheless,  to  obtain 
compatibility  and  in  preparation  for  performing  computations  that  are  more 
compounded,  the  technique  described  in  Eq.  (1-14)  is  utilized  to  Fourier  transform 
Kon(^j£y)  to  From  Eqs.  (1-14)  and  (13a),  and  after  an  appropriate 

normalization,  one  obtains 

v^„(x,®)  ~ 

R 

e„(a))  [k-i(1  +  R)"^]  ^  Sli^k^)  [G^^{kr,a))  exp{-i(x-Xa)K]]  ; 

r  =  o 

e^{(0)  =  [{llK^kl^)  {{nlakp„)^ ;  (-1)^/2  =  (_/)  _ 

The  normalization  adopted  in  Eq.  (16)  yields  for  a  hybrid  cylinder  in  the  absence  of  fluid 
loading,  at  the  position  of  the  external  ring  drive;  i.e.,  at  (Ix-x^  I/ft)  =  0,  a  unit 
amplitude.  In  Fig.  8a  the  normalized  modal  response  v„„(x,<o)  is  displayed  as  a  function 
of  (x/ft)  for  the  normalized  frequency  {(olm^)  =  0.175.  The  computations  in  this 
figure  are  performed  under  standard  parametric  values  for  a  hybrid  cylinder. 
Computations  pertaining  to  changes  in  the  fluid  loading  parameter  (e^)  from  the 
standard  value  of  10"^  to  10"^  and  10"^  are  superposed  on  the  original  Fig.  8a. 
[cf.  Figs.  2a-c.]  The  influence  of  the  fluid  surface  mass  loading  is  revealed  in  this  figure; 
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the  curves  pertaining  to  =  10  ^  and  10  ^  substantially  overlap,  whereas  the  curve 
for  =  10"^  lies  just  below,  but  consistently  adjacent.  There  is,  however,  no  evidence 
of  fluid  free  waves  in  Fig.  8a;  the  dominance  of  the  flexural  free  waves  appear 
overwhelming.  This  figure  is  repeated  in  Figs.  9a  and  10a,  except  that  a  natural  cylinder 
is  substituted  for  the  hybrid  cylinder  and  the  mode  index  (n)  is  selected  to  be  unity  and 
zero;  respectively.  Figures  9a  and  10a  are  identical  to  Fig.  8a  and,  therefore,  it  appears 
that  the  dominance  of  the  flexural  free  waves  overwhelms  not  only  the  fluid  free  waves, 
but  also  the  membrane  free  waves.  To  emphasize  this  identity.  Fig.  8a  is  jointly 
designated  Figs.  9a  and  10a.  Can  the  flexural  free  waves  be  rendered  weak  enough  so 
that  those  other  types  of  free  waves,  which  are  faintly,  but  clearly  discernible  in  Fig.  2 
and  3,  may  be  boldly  exposed  and  investigated?  Figures  2f  and  3f  and  g  suggest  that 
damping  the  flexural  free  waves  may  induce  the  desired  weakening  in  these  free  waves. 
Figures  8a,  9a  and  10a  are  repeated  in  Figs.  8b,  9b  and  10b,  respectively,  except  that  the 
flexural  loss  factor  {7]^)  is  changed  from  the  standard  value  of  10“^  to  3  x  10"^,  a 
substantial  increase  in  damping.  This  increase  in  damping  of  the  flexural  free  waves 
appears  more  dramatic  and  revealing  in  Figs.  8b,  9b  and  10b  than  in  the  corresponding 
Figs.  2d,  3d  and  3g  in  the  k  —domain.  In  Figs.  8b,  9b  and  10b  the  increase  spatial  decay 
of  the  flexural  free  waves  is  bottomed;  this  bottoming  is  not  due  to  lack  of 
computational  accuracy,  but,  rather,  due  to  the  existence  and  presence  of  other  types  of 
free  waves.  In  part,  to  help  substantiate  this  statement  Figs.  8c,  9c  and  10c  are  added.  In 
these  figures  the  fluid  loading  parameter  (e^)  is  set  equal  to  10"^  and,  in  conformity,  the 
flexural  loss  factor  (77^)  is  retained  equal  to  3  x  10"^.  Figures  8c,  9c  and  10c  are  also 
superposed  on  Figs.  8b,  9b  and  10b.  The  presence  of  fluid  free  waves  emerges  in  Fig.  8b 
and  the  presence  of  fluid  and  membrane  free  waves  emerge  in  Figs.  9b  and  10b. 
Moreover,  manifestation  of  interferences  among  the  various  types  of  free  waves  are 
discernible  in  these  figures.  Interferences  occur  in  "regions"  in  which  more  than  one 
type  of  free  waves  are  present  with  strengths  within  an  order  of  magnitude.  Thus,  in 
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Figs.  8b  and  c  interference  occurs  in  the  region  of  the  bottoming  knee;  this  interference 
is  between  the  flexural  and  fluid  free  waves.  In  this  region  the  strengths  of  these  two 
types  of  free  waves  are  within  an  order  of  magnitude.  Increasing  the  spatial  distance 
(}x-Xa\lb)  beyond  the  knee,  leaves  the  fluid  free  waves  alone;  the  flexural  free  waves 
are  weakened  out  by  the  high  damping.  Changes  in  the  strength  of  the  fluid  free  waves 
are  then  exhibited  in  Fig.  8b;  the  knees  appropriately  adjust  to  these  changes  over 
several  orders  of  magnitude.  In  passing  it  is  recalled  that  in  Fig.  8a  the  fluid  loading 
significantly  influences  the  surface  impedance  of  the  cylinder  only  when  >  10“^. 
Indeed,  this  influence  causes  the  saturation  and  reversal  in  the  strength  of  fluid  free 
waves  as  the  fluid  loading  parameter  (e^)  is  increased  through  10“^  on  to  10~^  The 
apparent  anomaly  is  associated  with  Le  Chartlier's  principle  [2].  In  Figs.  9b  and  c  and 
10b  and  c,  the  membrane  free  waves  join  in;  in  Figs.  9b  and  c  the  membrane  free  waves 
are  comprised  of  two  components;  longitudinal  and  shear  free  waves  and  in  Figs.  10b 
and  c  of  one  component;  longitudinal  free  waves  only.  These  differences  in  the  two  sets 
of  figures  help  decipher  the  free  waves  that  are  involved  in  the  interferences.  For 
example,  it  is  clear  from  Figs.  9c  and  10c  that  the  strengths  of  the  longitudinal  and  shear 
free  waves,  at  =  0.175,  are  within  an  order  of  magnitude  and,  therefore  they 

significantly  interfere  in  Fig.  9c,  whereas  this  interference  is  absent  in  Fig.  10c.  In  these 
two  figures  the  fluid  free  waves  do  not  participate  in  the  interference;  the  fluid  free 
waves  are  rendered  weak  by  selecting  the  fluid  loading  parameter  (e^)  equal  to  10"^  in 
these  figures.  Clearly  the  fluid  free  waves  participate  in  the  interferences  depicted  in 
Figs.  9b  and  10b  for  a  fluid  loading  parameter  >  lO'^;  in  particular  with  the 
membrane  free  waves  at  distances  Qx-x^l/b)  beyond  the  bottoming  knees.  The 
details  that  are  captured  in  the  interferences  in  these  figures  attest  to  the  accuracy 
attained  in  the  computations  via  Eq.  (1-14),  The  simplicity  of  the  summand  and, 
therefore,  the  simplicity  of  the  computations  performed  to  generate  Figs.  8  through  10, 
may  not  convincingly  constitute  a  critical  proof  of  the  viability  of  Eq.  (1-14). 
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Nonetheless,  the  kind  of  figures  just  presented  are  of  intrinsic  value,  especially  as  they 
relate  and  assist  in  deciphering  figures  displaying  more  compounded  computations; 
those  that  account  for  the  attachment  of  ribs  on  the  cylinder. 
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IV.  MODAL  DRIVE  IN  LIEU  OF  THE  RIBS  IN  THE  SPATIAL  DOMAIN 
The  normalized  modal  drive  P^{k,(o)  in  lieu  of  the  ribs  is  typically  depicted  in  the 
k  -domain  for  a  hybrid  and  a  natural  cylinder  in  Figs.  4  and  5,  respectively.  These 
figures  are  computed  under  standard  parametric  values  and  for  a  selected  variety  in 
these  values.  These  figures  demonstrate  that  this  quantity  is  aliased  in  {k)  with  respect 
to  (k’j);  the  expression  for  this  quantity  in  Eq.  (14)  confirms  this  aliasing.  Moreover,  the 
expression  for  P^^(k,co)  indicates  that  Eq.  (1-21)  suitability  befits  the  formula  for 
performing  the  Fourier  transformation  for  this  quantity.  Designating  p^{x,(o)  the 
Fourier  transform  of  P^(k,co),  Eqs.  (1-14),  (1-21)  -  (1-23)  and  (14)  are  used  to  derive 

=  Sh(x)  [5(x)  p^(x,co)}  ; 

Psn(x,(0)  =  [p,(x,CO)/ P^^((0)]  ,  (17) 

where 

P°n(x,0})  -r 

R 

(1  +  ^)"^  2  exp(-ixkr)  ll^{kr,(o)  S^ik,)  expiix^k,)}  ,  (18a) 

r  =  o 

lL(.k,<o)  =  [1  +  Sl(k)  Q^(k,co)r'  ,  (18b) 

and  it  is  recognized  that  in  Eqs.  (17)  and  (18)  there  exist  a  factorial  saving  of  (2/  + 1)^  in 
computational  steps  as  compared  with  estimating  p^(x,co)  from  P^^(k,(o)  employing 
the  discrete  Fourier  transform  procedure,  as  stated  in  Eq.  (l-2c).  The  absolute  values  of 
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p^n(x,co)  as  a  function  of  the  normalized  frequency  (o) I co^)  for  a  hybrid  and  a  natural 
cylinder  are  presented  in  Figs.  1 1  and  12,  respectively.  The  computations  are  performed 
for  (x/b)  =>  {Xj  lb)  =  7,14,  21  and  35  and  for  a  mode  index  (n)  equal  to  unity. 
Each  of  these  figures  is  composited  of  three  values  of  the  fluid  loading  parameter; 
namely,  -  10“^,  10“^  and  10“^.  The  presence  of  the  phenomenon  of  pass  and  stop 
bands  with  respect  to  the  flexural  free  waves  emerges  in  these  figures.  The  influence  of 
fluid  loading  on  this  phenomenon  is  also  distinguished  in  these  figures.  In  particular 
there  is  a  shift  to  lower  frequencies  of  the  stop  bands  with  increase  in  fluid  loading;  the 
shift  is  more  pronounced  the  higher  the  fluid  loading  is.  This  shift  in  frequencies  is 
related  to  the  added  surface  mass  component  in  the  fluid  surface  impedance  as  already 
discussed  with  respect  to  Figs.  2-5.  This  and  other  complemental  features  in  P^{k,(o) 
and  p^^(x,co)  are  facilitated  by  superimposing  the  corresponding  computations 
performed  on  each  quantity  in  the  same  figure.  Thus,  Figs.  4  and  5  simultaneously 
exhibit  both  quantities.  It  is  observed  in  these  figures  that  the  high  spatial  decay  in  the 
modal  drive  p^(x,0))  in  a  (flexural)  stop  band  in  the  x— domain,  is  appropriately 
commensurate  with  the  uniformity  in  the  drive  P^^(k,co)  in  a  stop  band  in  the 
^-domain.  This  uniformity  is  derived  from  cancellations  between  pairs  of  flexural 
dispersive  aliasing  orders  in  those  spectral  regions  in  which  they  cross.  In  the  absence 
of  these  cancellations,  as  is  the  case  in  a  first  order  model  in  which  interactions  among 
the  ribs  via  the  cylindrical  shell  are  neglected,  pass  and  stop  bands  are  also  absent  [2,13]. 
The  pass  and  stop  bands  are  largely  dominated  by  the  flexural  free  waves;  the  other 
types  of  free  waves  are  poorly  manifested  by  a  few  minor  kinks  in  the  curves  depicting 
p^n(x,(o)  and  a  few  corresponding  kinks  in  the  curves  depicting  P^^(k,(o).  Other 
features  in  these  quantities  can  be  similarly  related  from  the  k  -  domain  to  the  x  -  domain 
and  vice  versa;  e.g.,  the  shifts  in  frequencies  of  the  stop  bands  due  to  a  change  in  the 
character  of  the  ring  impedance  of  the  ribs  from  mass  to  stiffness  controlled  correctly 
track  in  the  complemental  representations  in  Figs.  4d  and  5d.  There  are,  of  course. 
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features  that  are  more  indigenous  to  one  domain  than  another.  The  aliasing  in  P^{k,(o) 
are  clearly  a  ^-domain  phenomenon.  On  the  other  hand,  in  a  stop  band,  when  the 
distance  (I  Xj  -x^\  lb)  exceeds  a  specific  value,  the  decay  in  p^^(x,co)  is  bottomed.  The 
bottomings  in  the  (flexural)  stop  bands  are  clearly  a  x- domain  phenomenon.  The 
bottoming  phenomenon  in  a  stop  band  is  reminiscent  of  the  same  named  phenomenon 
induced  by  the  increase  in  damping  of  the  flexural  free  waves  in  the  absence  of  ribs.  Is 
there  a  commonality  in  these  two  phenomena? 
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V.  RESPONSE  OF  A  REGULARLY  RIBBED  CYLINDER  IN  THE 

SPATIAL  DOMAIN 

The  normalized  modal  response  V^(k,co)  is  composed  of  two  terms,  V^{k,(o) 
and  as  stated  in  Eqs.  (12)  and  (13),  Formally,  from  Eqs.  (1-la),  (12)  and  (13) 

one  obtains 

v^(x,co)  =  v^^(x,Q))  -  v^(x,(o)  ,  (19) 

=  ^n(®)  (2;r)"^^^  I  dk  exp(-ixk)  V^^(k,co)  ,  (20a) 

V5„(^.to)  =  «n(to)  (2;r)"^^^  J  dk  exp(-ixk)  V^^(k,(o)  ,  (20b) 

where  e^(co)  is  the  normalization  factor  defined  in  Eq.  (16),  The  Fourier  transform  of  the 
first  term  is  dealt  with  in  Section  HI;  in  this  section  the  second  term  is  added  to  the 
computational  task.  As  already  discussed  in  the  A: -domain,  neither  V^{k,(o)  nor 
Vg„(k,co)  are  strictly  aliased.  In  the  first  term  neither  of  its  two  factors  are  aliased  and, 
therefore,  no  factorial  saving  can  be  accrued  from  the  form  of  its  evaluation.  In  the 
second  term  the  factor  G„^{k,(o)  is  not  aliased;  however,  the  second  factor,  being  the 
modal  drive  P^(k,co)  in  lieu  of  the  ribs,  is  aliased.  Identifying  G{k,(o)  as  G^^(k,a})  and 
K(k,a))  as  P^^(k,co)  in  Eq,  (1-14),  and  with  Eqs,  (l-16b)  and  (14a)  in  tow,  a  factorial 
saving  of  (2J  -f-1),  and  with  Eqs,  (l-19b)  and  (14b)  in  tow,  a  factorial  saving  of  (2/-f  1)^ 
of  computational  steps  may  be  accrued  in  the  evaluation  of  v^„(x,  (o)  in  Eq.  (20b).  [This 
kind  of  saving  is  not  afforded  using  Eq.  (l-2c).]  From  Eqs.  (1-14),  (l-19b),  (14b)  and 
(20b)  one  obtains 
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e„(co)  [k:i(1  + /?r^]  li{.K,(o)  S^ik,)  [Q^nik^co)  exp(ix^k^)} 

r  =  o 

Sl(k,)  [exp(-ixk,)  G„(k,,m)]  ; 

e„(o)  =  l(.2/z^k%)  Unlak^^f-l]]!^'^  ■  (-!)'«  =  (-0  .  (21) 

The  normalization  and  evaluation  of  the  second  term  in  Eq.  (19)  are  rendered  compatible, 
in  Eq.  (21),  with  those  of  the  first  term  in  Eq.  (16).  A  question  may  be  raised:  Since  the 
aliasings  phenomenon  in  P^(k,co)  is  spoiled  in  \^(^,<»)  and  the  pass  and  stop  bands 
phenomenon  is  revealed  in  the  modal  drive  p^^(x,co),  is  the  pass  and  stop  bands 
phenomenon  spoiled  in  the  modal  response  v„(jc,to)?  In  other  words,  how  akin  are 
these  two  phenomena?  The  absolute  value  of  the  normalized  modal  response  v^(x,(o)  is 
displayed  as  a  function  of  (x/b),  for  a  hybrid  cylinder  in  Fig.  13  and  a  natural  cylinder  in 
Figs.  14  and  15;  in  Figs.  14  and  15  the  mode  index  is  set  equal  to  unity  and  zero, 
respectively.  In  Figs.  13a,  14a  and  15a  the  frequency  is  set  at  (col(0^)  =  0.175  with 
the  fluid  loading  parameter  =  10”^,  lO""^,  10“^  and  10"^  and  at  (colco^)  =  0.158 
with  =  10“^  It  is  observed  that  data  pertaining  to  the  frequency  (to  /  )  =  0.175 

with  =  lO"^,  10”^,  and  I0~^  and  the  frequency  (co/to^)  =  0.158  with  =  10~^ 
lie  in  a  (flexural)  stop  band.  The  data  pertaining  to  a  frequency  (co/co^)  =  0.175  with 
Cc  =  10"^  substantially  lie  in  a  pass  band;  the  shift  in  the  frequencies  of  the  stop  bands 
under  the  influence  of  heavy  fluid  loading  accounts  for  this  feature  in  the  data  with 
respect  to  =  10”^  [cf.Figs.il  and  12.]  Figures  13a,  14a  and  15a  convincingly 
indicate  that  the  pass  and  stop  bands  phenomenon  is  not  closely  akin  to  the 
phenomenon  of  aliasing;  the  spoiling  of  the  aliasing  in  V^(k,(o)  does  not  entail  the 
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spoiling  of  the  pass  and  stop  bands  in  v„(jc,co).  It  emerges,  however,  that  data  in 
Figs.  13a,  14a  and  15a  that  lie  in  the  stop  band  at  {col (O^)  =  0.175  with  -  10“^, 
10"^  and  10“^  and  at  (colco^)  =  0.158  with  =  10"\  are  reminiscent  of  the 
corresponding  data  in  Figs.  8b,  9b  and  10b,  respectively.  This  reminiscence  holds  even 
with  respect  to  the  reversal  of  strength  between  the  fluid  free  waves  with  =  10“^ 
and  10"^;  the  anomaly  is  particularly  discernible  in  Figs.  8b  and  13a.  In  the  vein  of 
Figs.  8c,  9c  and  10c,  the  curves  in  Figs.  13a,  14a  and  15a  that  pertain  to  a  fluid  loading 
parameter  =  10“^  are,  in  isolation,  displayed  in  Figs.  13b,  14b  and  15b,  respectively. 
Notwithstanding  that  there  exist  minor  and  accountable  differences  in  levels,  the 
resemblance  between  Figs.  13a  and  b,  14a  and  b,  and  15a  and  b  and  Figs.  8b  and  c,  9b 
and  c,  and  10b  and  c,  respectively,  is  striking.  The  resemblance  even  matches  the 
bottoming  and  interference  phenomena  described  earlier.  Indeed,  comparison  among 
the  correspond  figures  in  the  two  sets  leads  one  to  conclude  that  the  regularly  attached 
ribs  in  a  stop  band  play  a  role  that  is  commensurate  with  the  role  played  by  the  increase 
in  flexural  damping  in  the  absence  of  ribs.  This  conclusion  does  not  take  account  of  the 
undulations  in  level  within  a  bay  that  exist  in  the  former  set,  but  are  absent  in  the  latter 
set  of  figures  [14].  The  undulations  are  avoided  in  the  former  set  of  figures  by  an 
appropriate  selection  of  the  position  (jc^  /  b)  of  the  external  drive  and  the  regularity  of 
the  sampling  positions  (x/b)  along  the  axis  of  the  cylinder  [14].  These  undulations  are 
not  significant  enough,  however,  to  altercate  the  striking  resemblance  between  these 
two  sets  of  figures. 

Figures  13a  and  b,  14a  and  b,  and  15a  and  b  are  made  possible  only  with  an 
accuracy  and  a  reproductive  details  that  transcend  several  orders  of  magnitude.  In  the 
computations  that  underlie  these  figures  the  overall  wavenumber  span  (k,^)  and  the 
sampling  wavenumber  span  (kTi)  are  related  largely  by  a  7  =  50,  yielding  a  factorial 

0  A 

saving  of  10  and,  in  some  cases,  of  10  of  computational  steps  in  computations 
performed  with  the  double-sum  technique  versus  computations  performed  with 
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techniques  that  are  directly  akin  to  the  discrete  Fourier  transform  procedures.  Whether 
this  kind  of  factorial  savings  is  cmcial  to  the  generation  of  the  figures  presented  in  this 
report  is  a  matter  of  memory  and  patience.  The  former  can  be  bought  and  the  latter  can 
be  acquired  by  practice,  notwithstanding  that  shortages  on  both  accounts  are  thinning 
out.  Nonetheless,  that  the  double-sum  technique  has  been  successfully  employed  to 
furnish  the  figures  just  presented  is  by  now  a  fact. 


Fluid  (p,  c)  (r,  0,  ([))  (r,  x,  (b) 
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Fig.  1.  A  sketch  of  a  regularly  ribbed  cylinder  and  coordinate  system. 
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Normalized  modal  response  of 


dejceutd 


d^jwgutd 


VnBar2p 


Fig.  6.  Normalized  response  V„  {k,  ft))  of  a 
regularly  ribbed  hybrid  cylinder  as  a 
function  of  the  normalized  axial 
wavenumber  {ak)  in  a  frequency 
waterfall  format. 

a.  Under  standard  values,  [cf.  Eq.  (15).] 

a.  1 .  Clipping  of  Fig.  6.a. 

b.  =  10“^. 

c.  Be  =  10“^ 

d.  Ribs’  ring  impedance  stiffness  controlled. 

e.  n  =  0. 
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36  42  4S5460667278S4 


-18-12-6  0  6  12  18  24  30  36  42  48  54  60  66  72  78  84 

(xtb) 


b) 

5a  Fig.  15b 

5  changed  from  the  standard  value  of  unity  to  zero. 

1  =  10"^,  lO-"*,  10“^  and  10“^  and  (co/o)^)  =  0.158  with  =  10"^ 

[t  =  10“^. 
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